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Using a data sample corresponding to an integrated luminosity of 0.35 fb _1 collected by LHCb 
in 2011, we report the first evidence of CP violation in the decays of B® mesons to K ± n^ 1 pairs, 
A C p{B ( l -> Kir) = 0.27 ± 0.08 (stat) ± 0.02 (syst), with a significance of 3.3a. Furthermore, we 
report the first observation of CP violation in B° decays at a hadron collider, Acp(B° — > Ktv) — 
—0.088 ± 0.011 (stat) ± 0.008 (syst), with a significance exceeding 6a. 
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The violation of CP symmetry, i.e. the non-invariance 
of fundamental forces under the combined action of the 
charge conjugation (C) and parity (P) transformations, 
is well established in the K° and B° meson systems pOh 
H] . Recent results from the LHCb collaboration have also 
provided evidence for CP violation in the decays of D° 
mesons [5 . Consequently, there now remains only one 
neutral heavy meson system, the B®, where CP viola- 
tion has not yet been seen. All current experimental 
measurements of CP violation in the quark flavor sector 
are well described by the Cabibbo-Kobayashi-Maskawa 
mechanism [5] [7] which is embedded in the framework of 
the Standard Model (SM). However, it is believed that 
the size of CP violation in the SM is not sufficient to ac- 
count for the asymmetry between matter and antimatter 
in the Universe [8], hence additional sources of CP sym- 
metry breaking are being searched for as manifestations 
of physics beyond the SM. 

In this Letter we report measurements of direct CP 
violating asymmetries in B° — > K + tt~ and B° s — > K~ir + 
decays using data collected with the LHCb detector. The 
inclusion of charge-conjugate modes is implied except in 
the asymmetry definitions. CP violation in charmless 
two-body B decays could potentially reveal the presence 
of physics beyond the SM [9TH3"]. and has been exten- 
sively studied at the B factories and at the Tevatron [14T - 
[16] . The direct CP asymmetry in the B (B®) decay rate 
to the final state / = K + n~ (f s = K~ix + ) is defined as 



.4 



CP 



4> 



L B 



(<0 



/(s 



B (s) 



f(s) 



(1) 



where $[A, Y] = (X — Y)/{X + Y) and / (s) denotes the 
charge-conjugate of ft 8 y 

LHCb is a forward spectrometer covering the pseudo- 
rapidity range 2 < r\ < 5, designed to perform flavor 
physics measurements at the LHC. A detailed description 
of the detector can be found in Ref. [17 . The analysis 
is based on pp collision data collected in the first half of 
2011 at a center-of-mass energy of 7 TeV, corresponding 
to an integrated luminosity of 0.35 fb _1 . The polarity of 
the LHCb magnetic field is reversed from time to time in 
order to partially cancel the effects of instrumental charge 
asymmetries, and about 0.15 fb _1 were acquired with one 
polarity and 0.20 fb _1 with the opposite polarity. 

The LHCb trigger system comprises a hardware trig- 
ger followed by a High Level Trigger (HLT) implemented 
in software. The hadronic hardware trigger selects high 
transverse energy clusters in the hadronic calorimeter. A 
transverse energy threshold of 3.5 GeV has been adopted 
for the data set under study. The HLT first selects events 
with at least one large transverse momentum track char- 
acterized by a large impact parameter, and then uses al- 
gorithms to reconstruct D and B meson decays. Most of 
the events containing the decays under study have been 
acquired by means of a dedicated two-body HLT selec- 
tion. To discriminate between signal and background 



TABLE I. Summary of selection criteria adopted for the mea- 
surement of A C p(B° -> Kty) and A C p{B° -> Kir). 



Variable 


A C p(B° -> Kty) 


Acp(B° -+ 


Track quality xV n df 


< 3 


< 3 


Track p T [GeV/c] 


> 1.1 


> 1.2 


Track dip [mm] 


> 0.15 


> 0.20 


max(pf, p£) [GeV/c] 


> 2.8 


> 3.0 


max(df P , dip) [mm] 


> 0.3 


> 0.4 


c?ca [mm] 


< 0.08 


< 0.08 


pf [GeV/c] 
dfp [mm] 


> 2.2 


> 2.4 


< 0.06 


< 0.06 


tnn [ps] 


> 0.9 


> 1.5 



events, this trigger selection imposes requirements on: 
the quality of the online-reconstructed tracks (\ 2 per 
degree of freedom), their transverse momenta (pt) an d 
their impact parameters (dip, defined as the distance be- 
tween the reconstructed trajectory of the track and the 
pp collision vertex); the distance of closest approach of 
the daughter tracks (c2ca); the transverse momentum of 
the B meson candidate fj*j ), its impact parameter (df P ) 
and its decay time (t vv , calculated assuming the decay 
into 7r + 7r~). Only B candidates within the ttit invariant 
mass range 4.7-5.9 GeV/c 2 are accepted. The tttt mass 
hypothesis is conventionally chosen to select all charmless 
two-body B decays using the same criteria. 

Offline selection requirements are subsequently ap- 
plied. Two sets of criteria have been optimized with 
the aim of minimizing the expected uncertainty either on 
A C p(B° -> A"7r) or on A C p(B° -> Kit). In addition to 
more selective requirements on the kinematic variables 
already used in the HLT, two further requirements on 
the larger of the transverse momenta and of the impact 
parameters of the daughter tracks are applied. A sum- 
mary of the two distinct sets of selection criteria is re- 
ported in Table [i] In the case of B® —> Kir decays a 
tighter selection is needed because the probability for a 
b quark to decay as B® — > Kty is about 14 times smaller 
than that to decay as B° — » A7r [T5], and consequently 
a stronger rejection of combinatorial background is re- 
quired. The two samples passing the event selection are 
then subdivided into different final states using the parti- 
cle identification (PID) provided by the two ring-imaging 
Cherenkov (RICH) detectors. Again two sets of PID se- 
lection criteria are applied: a loose set optimized for the 
measurement of Acp(B q —> Kir) and a tight set for that 
of A CP (B° -> Ktt). 

To estimate the background from other two-body B 
decays with a misidentified pion or kaon (cross-feed back- 
ground), the relative efficiencies of the RICH PID selec- 
tion criteria must be determined. The high production 
rate of charged D* mesons at the LHC and the kine- 



matic characteristics of the D* 



decay 



chain make such events an ideal calibration sample for 
the PID of kaons and pions. In addition, for calibrating 
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KV invariant mass (GeV/c 2 ) K jt + invariant mass (GeV/c 2 ) 



FIG. 1. Invariant Kit mass spectra obtained using the event selection adopted for the best sensitivity on (a, b) Acp(B° — > Kir) 
and (c, d) Acp{B ( l — > Kn). Plots (a) and (c) represent the K + ir~ invariant mass whereas plots (b) and (d) represent the 
K~n + invariant mass. The results of the unbinned maximum likelihood fits are overlaid. The main components contributing 
to the fit model are also shown. 



the response of the RICH system for protons, a sample of 
A pir~ decays is used. PID information is not used to 
select either sample, as the selection of pure final states 
can be realized by means of kinematic criteria alone. The 
production and decay kinematics of the D° —> K~w + 
and A — > pir~ channels differ from those of the B de- 
cays under study. Since the RICH PID information is 
momentum dependent, the distributions obtained from 
calibration samples are reweighted according to the mo- 
mentum distributions of B daughter tracks observed in 
data. 

Unbinned maximum likelihood fits to the Kit mass 
spectra of the selected events are performed. The B° — » 
Kit and B® — > Kir signal components are described 
by single Gaussian functions convolved with a function 
which describes the effect of final state radiation on the 
mass lineshape [19] . The background due to partially 
reconstructed three-body B decays is parameterized by 
means of an ARGUS function [5D] convolved with a Gaus- 
sian resolution function. The combinatorial background 
is modeled by an exponential and the shapes of the 
cross-feed backgrounds, mainly due to B — > ir + ir~ and 
B° -> K+K~ decays with one misidentified particle in 
the final state, are obtained from Monte Carlo simula- 
tions. The B° -> ir+ir- and B° -> K + K~ cross-feed 



background yields are determined from fits to the ir + ir~ 
and K + K~ mass spectra respectively, using events se- 
lected by the same offline selection as the signal and tak- 
ing into account the appropriate PID efficiency factors. 
The K + ir~ and K~ir + mass spectra for the events pass- 
ing the two offline selections are shown in Fig. [I] 

From the two mass fits we determine respectively 
the signal yields N(B° -> Kit) = 13250 ± 150 and 
N(B° -> Ktt) = 314 ± 27, as well as the raw yield 
asymmetries A law (B° -> Ktt) = -0.095 ± 0.011 and 
A raw (£° -> Kit) = 0.28 ± 0.08, where the uncertain- 
ties are statistical only. In order to determine the CP 
asymmetries from the observed raw asymmetries, effects 
induced by the detector acceptance and event reconstruc- 
tion, as well as due to strong interactions of final state 
particles with the detector material, need to be taken 
into account. Furthermore, the possible presence of a 
B9s — B( s j production asymmetry must also be consid- 
ered. The CP asymmetry is related to the raw asymme- 
try by Acp — A law — A a, where the correction A a is 
defined as 

A A (B° {S) -> Ktt) = Cd(s)A B (K7r) + K d(s) A P (B° [s) ), (2) 

where Cd — 1 an( i Cs — ~ U following the sign con- 
vention for / and f s in Eq. ([T]). The instrumental 
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asymmetry Ad(Ktt) is given in terms of the detection 
efficiencies £d of the charge-conjugate final states by 
A D (Kn) = $[e D (i^~7r + ), £d(A+7t~)], and the produc- 
tion asymmetry Ap(B? a \) is defined in terms of the 

production rates, R(B^) and R(B^), as 



* W and 

A P (B°s) = $[i?(B" s) ), R(B9 s) )]. The factor n d(s) takes 
into account dilution due to neutral B^ meson mixing, 
and is defined as 



Jo 



t cos(Am d{s) t)e(Bl s) ^KTr; t)dt 



/ oo e-^(^cosh(^>t) £ (i?0 s) - 



■Ktt; t)dt 



(3) 



where e(B° — > A-7r; t) and e(-Bg — >• Kit; t) are the accep- 
tances as functions of the decay time for the two recon- 
structed decays. To calculate n d and n s we assume that 
Ar^ = and we use the world averages for T^, Amj, T s , 
Am s and Ar s [3]. The shapes of the acceptance func- 
tions are parameterized using signal decay time distribu- 
tions extracted from data. We obtain K d = 0.303 ±0.005 
and k s = —0.033 ± 0.003, where the uncertainties are 
statistical only. In contrast to n d , the factor n s is small, 
owing to the large oscillation frequency, thus leading 
to a negligible impact of a possible production asymme- 
try of Bg mesons on the corresponding CP asymmetry 
measurement. 

The instrumental charge asymmetry Ay>(Kit) can 
be expressed in terms of two distinct contributions 
A-d(Kit) = Ai(Ktt) + a(Kir)An,(Kir), where Ai(K%) 
is an asymmetry due to the different strong interaction 
cross-sections with the detector material of K + ir~ and 
K~ir + final state particles, and An(Kn) arises from the 
possible presence of a reconstruction or detection asym- 
metry. The quantity Ai(Kit) does not change its value by 
reversing the magnetic field, as the difference in the inter- 
action lengths seen by the positive and negative particles 
for opposite polarities is small. By contrast, Ar(Kit) 
changes its sign when the magnetic field polarity is re- 
versed. The factor a(Kir) accounts for different signal 
yields in the data sets with opposite polarities, due to 
the different values of the corresponding integrated lumi- 
nosities and to changing trigger conditions in the course 
of the run. It is estimated by using the yields of the 
largest decay mode, i.e. B° — » Kit, determined from the 
mass fits applied to the two data sets separately. We ob- 
tain a{Kn) = <5>[N u p(B° -> Kit), N down (B° -> Kit)] = 
-0.202 ± 0.011 , where "up" and "down" denote the di- 
rection of the main component of the dipole field. 

The instrumental asymmetries for the final state Kir 
are measured from data using large samples of tagged 
D*+ -> D a (K-TT+)ir+ and D*+ -> D°(K-K+)ir+ de- 
cays, and untagged D° —> K~n + decays. The combina- 
tion of the integrated raw asymmetries of all these decay 
modes is necessary to disentangle the various contribu- 
tions to the raw asymmetries of each mode, notably in- 
cluding the Kir instrumental asymmetry as well as that of 
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FIG. 2. Distributions of the invariant mass or invari- 
ant mass difference of (a) D° — > K~n + , (b) D* + — > 
D°(K-tv + )it + , (c) D*+ -» D°(K-K + )ty+ and (d) B° -> 
J/i/j(h + fi~)K*° (K + -k~). The results of the maximum like- 
lihood fits are overlaid. 



the pion from the D* + decay, and the production asym- 
metries of the D* + and D° mesons. In order to determine 
the raw asymmetry of the D° — » Kit decay, a maximum 
likelihood fit to the K~ir + and K + ir~ mass spectra is 
performed. For the decays D* + — > D a (K~-K + )n + and 
D* + — > D°(K~K + )ir + , we perform maximum likelihood 
fits to the discriminating variable 8m = Md* — Moo, 
where Mo* and Mjyo are the reconstructed D* and D° 
invariant masses respectively. Approximately 54 million 

+ W+ 



D° 



decays, 7.5 million D*+ D°(A^7r+)7r 



and f.l million D* + -» D°(K~ K + )ir + decays are used. 
The mass distributions are shown in Fig. [2] (a), (b) 
and (c). The D° — > K~ir + signal component is mod- 
eled as the sum of two Gaussian functions with com- 
mon mean convolved with a function accounting for fi- 
nal state radiation 19J , on top of an exponential com- 
binatorial background. The D* + D°(K^ir + )ir + and 
D* + —> D°(K~ K + )ir + signal components are modeled 
as the sum of two Gaussian functions convolved with a 
function taking account of the asymmetric shape of the 
measured distribution [5]. The background is described 
by an empirical function of the form I — e -( Sm - 8m o)/i ; 
where 5mo and £ are free parameters. Using the cur- 
rent world average of the integrated CP asymmetry for 
the D° — > K~K + decay [21] and neglecting CP viola- 
tion in the Cabibbo-favored D° —> K~ir + decay |22) . 
from the raw yield asymmetries returned by the mass 
fits we determine Ai(Kir) = (-1.0 ± 0.2) x 10~ 2 and 
A^{Kir) = (—1.8 ± 0.2) x 10~ 3 , where the uncertainties 
are statistical only. 

The possible existence of a B° — B° production 
asymmetry is studied by reconstructing a sample of 
B° — s- J/ipK*° decays. CP violation in 6 ^ ccs 
transitions, which is predicted in the SM to be at 
the 10~ 3 level 23) . is neglected. The raw asymmetry 
A law (B° — > J/ipK*°) is determined from an unbinned 
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maximum likelihood fit to the J/ip(fJ- + Li )K*°(K + ir ) 
and J/ip(fJ- + /J-~)K (K~ir + ) mass spectra. The signal 
mass peak is modeled as the sum of two Gaussian func- 
tions with common mean, whereas the combinatorial 
background is modeled by an exponential. The data 
sample contains approximately 25 400 B° —> J/i/jK*° de- 
cays. The mass distribution is shown in Fig. [2] (d). To 
determine the production asymmetry we need to correct 
for the presence of instrumental asymmetries. Once the 
necessary corrections are applied, we obtain a value for 
the B° production asymmetry A P (B°) = 0.010 ± 0.013, 
where the uncertainty is statistical only. 

By using the instrumental and production asym- 
metries, the correction factor to the raw asymmetry 
A A (B° -> Kir) = -0.007 ± 0.006 is obtained. Since 
the B® meson has no valence quarks in common with 
those of the incident protons, its production asymmetry 
is expected to be smaller than for the B°, an expectation 
that is supported by hadronization models as discussed 
in Ref . |24j . Even conservatively assuming a value of the 
production asymmetry equal to that for the B°, owing 
to the small value of k s the effect of Ap(B®) is negligible, 
and we find A A (B° -> Kir) = 0.010 ± 0.002. 

The systematic uncertainties on the asymmetries fall 
into the following main categories, related to: (a) PID 
calibration; (b) modeling of the signal and background 
components in the maximum likelihood fits; and (c) in- 
strumental and production asymmetries. Knowledge of 
PID efficiencies is necessary in this analysis to compute 
the number of cross-feed background events affecting the 
mass fit of the B a —> Kir and B® — > Kir decay chan- 
nels. In order to estimate the impact of imperfect PID 
calibration, we perform unbinned maximum likelihood 
fits after having altered the number of cross-feed back- 
ground events present in the relevant mass spectra ac- 
cording to the systematic uncertainties affecting the PID 
efficiencies. An estimate of the uncertainty due to pos- 
sible imperfections in the description of the final state 
radiation is determined by varying, over a wide range, 
the amount of emitted radiation [19] in the signal line- 
shape parameterization. The possibility of an incorrect 
description of the core distribution in the signal mass 
model is investigated by replacing the single Gaussian 
with the sum of two Gaussian functions with a common 
mean. The impact of additional three-body B decays in 
the Kir spectrum, not accounted for in the baseline fit 
- namely B — > iririr where one pion is missed in the 
reconstruction and another is misidentified as a kaon — 
is investigated. The mass lineshape of this background 
component is determined from Monte Carlo simulations, 
and then the fit is repeated after having modified the 
baseline parameterization accordingly. For the modeling 
of the combinatorial background component, the fit is 
repeated using a first-order polynomial. Finally, for the 
case of the cross- feed backgrounds, two distinct system- 



TABLE II. Summary of systematic uncertainties on 
A C p(B° -> Kir) and A C p(B° -> Kir). The categories (a), 
(b) and (c) defined in the text are also indicated. The total 
systematic uncertainties given in the last row are obtained by 
summing the individual contributions in quadrature. 



Systematic uncertainty 


A C p(B° -> 


Kir) A C p(B° -> Kir) 


(a) PID calibration 


0.0012 


0.001 


' b *Final state radiation 


0.0026 


0.010 


' b ' Signal model 


0.0004 


0.005 


' b ' Combinatorial background 


0.0001 


0.009 


3-body background 


0.0009 


0.007 


' b ' Cross-feed background 


0.0011 


0.008 


' c 'lnstr. and prod. asym. (Aa) 


0.0078 


0.005 


Total 


0.0084 


0.019 



atic uncertainties are estimated: one due to a relative 
bias in the mass scale of the simulated distributions with 
respect to the signal distributions in data, and another 
accounting for the difference in mass resolution between 
simulation and data. All the shifts from the relevant 
baseline values are accounted for as systematic uncer- 
tainties. Differences in the kinematic properties of B de- 
cays with respect to the charm control samples, as well 
as different triggers and offline selections, are taken into 
account by introducing a systematic uncertainty on the 
values of the A A corrections. This uncertainty dominates 
the total systematic uncertainty related to the instru- 
mental and production asymmetries, and can be reduced 
in future measurements with a better understanding of 
the dependence of such asymmetries on the kinematics 
of selected signal and control samples. The systematic 
uncertainties for Acp(B° -> Kir) and Acp{B° s -» Kir) 
are summarized in Table [TTJ. 

In conclusion we obtain the following measurements of 
the CP asymmetries: 

A CP (B Q -> Kir) = -0.088 ± 0.011 (stat) ± 0.008 (syst) 

and 

A C p(B° s -> Kir) = 0.27 ± 0.08 (stat) ± 0.02 (syst). 

The result for Acp(B° — »• Kir) constitutes the most 
precise measurement available to date. It is in good 
agreement with the current world average provided 
by the Heavy Flavor Averaging Group Acp(B° — > 
Kir) = -0.098l^u EQ- Dividing the central value of 
Acp(B q — > Kir) by the sum in quadrature of the sta- 
tistical and systematic uncertainties, the significance of 
the measured deviation from zero exceeds 6a, making 
this the first observation (greater than 5cr) of CP vio- 
lation in the B meson sector at a hadron collider. The 
same significance computed for Acp(B® — > Kir) is 3.3c, 
therefore this is the first evidence for CP violation in the 
decays of B® mesons. The result for Acp(B® — > Kir) 
is in agreement with the only measurement previously 
available [Tfj] , 
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